A Genuinely Hybrid, Multiscale 3D Cancer Invasion and Metastasis Modelling Framework.

We introduce in this paper substantial enhancements to a previously proposed hybrid multiscale cancer invasion modelling framework to better reflect the biological reality and dynamics of cancer. These model updates contribute to a more accurate representation of cancer dynamics, they provide deeper insights and enhance our predictive capabilities. Key updates include the integration of porous medium-like diffusion for the evolution of Epithelial-like Cancer Cells and other essential cellular constituents of the system, more realistic modelling of Epithelial-Mesenchymal Transition and Mesenchymal-Epithelial Transition models with the inclusion of Transforming Growth Factor beta within the tumour microenvironment, and the introduction of Compound Poisson Process in the Stochastic Differential Equations that describe the migration behaviour of the Mesenchymal-like Cancer Cells. Another innovative feature of the model is its extension into a multi-organ metastatic framework. This framework connects various organs through a circulatory network, enabling the study of how cancer cells spread to secondary sites.

Stochastic differential equation modelling of cancer cell migration and tissue invasion.

Invasion of the surrounding tissue is a key aspect of cancer growth and spread involving a coordinated effort between cell migration and matrix degradation, and has been the subject of mathematical modelling for almost 30 years. In this current paper we address a long-standing question in the field of cancer cell migration modelling. Namely, identify the migratory pattern and spread of individual cancer cells, or small clusters of cancer cells, when the macroscopic evolution of the cancer cell colony is dictated by a specific partial differential equation (PDE). We show that the usual heuristic understanding of the diffusion and advection terms of the PDE being one-to-one responsible for the random and biased motion of the solitary cancer cells, respectively, is not precise. On the contrary, we show that the drift term of the correct stochastic differential equation scheme that dictates the individual cancer cell migration, should account also for the divergence of the diffusion of the PDE. We support our claims with a number of numerical experiments and computational simulations.

A novel 3D atomistic-continuum cancer invasion model: In silico simulations of an in vitro organotypic invasion assay.

We develop a three-dimensional genuinely hybrid atomistic-continuum model that describes the invasive growth dynamics of individual cancer cells in tissue. The framework explicitly accounts for phenotypic variation by distinguishing between cancer cells of an epithelial-like and a mesenchymal-like phenotype. It also describes mutations between these cell phenotypes in the form of epithelial-mesenchymal transition (EMT) and its reverse process mesenchymal-epithelial transition (MET). The proposed model consists of a hybrid system of partial and stochastic differential equations that describe the evolution of epithelial-like and mesenchymal-like cancer cells, respectively, under the consideration of matrix-degrading enzyme concentrations and the extracellular matrix density. With the help of inverse parameter estimation and a sensitivity analysis, this three-dimensional model is then calibrated to an in vitro organotypic invasion assay experiment of oral squamous cell carcinoma cells.

Energy intake functions and energy budgets of ectotherms and endotherms derived from their ontogenetic growth in body mass and timing of sexual maturation.

Ectothermic and endothermic vertebrates differ not only in their source of body temperature (environment vs. metabolism), but also in growth patterns, in timing of sexual maturation within life, and energy intake functions. Here, we present a mathematical model applicable to ectothermic and endothermic vertebrates. It is designed to test whether differences in the timing of sexual maturation within an animal's life (age at which sexual maturity is reached vs. longevity) together with its ontogenetic gain in body mass (growth curve) can predict the energy intake throughout the animal's life (food intake curve) and can explain differences in energy partitioning (between growth, reproduction, heat production and maintenance, with the latter subsuming any other additional task requiring energy) between ectothermic and endothermic vertebrates. With our model we calculated from the growth curves and ages at which species reached sexual maturity energy intake functions and energy partitioning for five ectothermic and seven endothermic vertebrate species. We show that our model produces energy intake patterns and distributions as observed in ectothermic and endothermic species. Our results comply consistently with some empirical studies that in endothermic species, like birds and mammals, energy is used for heat production instead of growth, and with a hypothesis on the evolution of endothermy in amniotes published by us before. Our model offers an explanation on known differences in absolute energy intake between ectothermic fish and reptiles and endothermic birds and mammals. From a mathematical perspective, the model comes in two equivalent formulations, a differential and an integral one. It is derived from a discrete level approach, and it is shown to be well-posed and to attain a unique solution for (almost) every parameter set. Numerically, the integral formulation of the model is considered as an inverse problem with unknown parameters that are estimated using a series of empirical data.

A Multiscale Approach to the Migration of Cancer Stem Cells: Mathematical Modelling and Simulations.

We propose a multiscale model for the invasion of the extracellular matrix by two types of cancer cells, the differentiated cancer cells and the cancer stem cells. We investigate the epithelial mesenchymal-like transition between them being driven primarily by the epidermal growth factors. We moreover take into account the transdifferentiation program of the cancer stem cells towards the cancer-associated fibroblast cells as well as the fibroblast-driven remodelling of the extracellular matrix. The proposed haptotaxis model combines the macroscopic phenomenon of the invasion of the extracellular matrix by both types of cancer cells with the microscopic dynamics of the epidermal growth factors. We analyse our model in a component-wise manner and compare our findings with the literature. We investigate pathological situations regarding the epidermal growth factors and accordingly propose "mathematical-treatment" scenarios to control the aggressiveness of the tumour.

An extended Filament Based Lamellipodium Model produces various moving cell shapes in the presence of chemotactic signals.

The Filament Based Lamellipodium Model (FBLM) is a two-phase two-dimensional continuum model, describing the dynamics of two interacting families of locally parallel actin filaments (Oelz and Schmeiser, 2010b). It contains accounts of the filaments' bending stiffness, of adhesion to the substrate, and of cross-links connecting the two families. An extension of the model is presented with contributions from nucleation of filaments by branching, from capping, from contraction by actin-myosin interaction, and from a pressure-like repulsion between parallel filaments due to Coulomb interaction. The effect of a chemoattractant is described by a simple signal transduction model influencing the polymerization speed. Simulations with the extended model show its potential for describing various moving cell shapes, depending on the signal transduction procedure, and for predicting transients between non-moving and moving states as well as changes of direction.